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ABSTRACT 

Recent calculations of coUisionless, Case-B, He I emissivities were performed 
by Bauman et al. (2005). The source code used in the calculation has been freely 
available online since that paper was published. A number of changes have been 
made to simplify the use of the code by third parties. Here I provide details on 
how to obtain, compile, and execute the program and interpret the results. 



1. Introduction 

This document serves a simple purpose: to inform interested workers about how to 
obtain and use the "J-resolved" computer program discussed in Bauman et al. (2005), 
hereafter referred to as "the code." There are no new scientific results presented here, and 
the methods and predictions of the code are unchanged. 

2. Obtaining the Code 

The code can be obtained at 'http:/ /www.pa.uky.edu/~rporter/j-resolved/". There are 
two versions available, but the difference is only in the packaging. The file heliumRecomb . tar . 
is convenient for use on Unix and Linux systems, while the file heliumRecomb.zip is more 
convenient for Microsoft Windows machines. 



3. Compiling the Code 

While I have tested the code on several systems, I cannot guarantee that it will compile 
correctly on any given system. Nonetheless, the code will likely work on any version of 
Windows, Unix, Linux, or BSD. I have made no attempt to compile the code in MacOS. 
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Before compiling, you may wish to edit the file Helium. h so that the variable 
OUTPUT_DIR on the first line points to a directory other than where the source files 
are located. Now you are ready to compile. The preferred compiler is the open-source 
GNU compiler, gcc (available at http://gcc.gnu.org/). Compile with gcc using the following 
commands: 



gcc -c *.c 

gcc -o EXECUTABLE *.o -Im 



where EXECUTABLE is whatever filename you wish to give to your executable 
(i.e., helium.exe). 

The code can also be compiled with Microsoft Visual Studio. First, create a new project 
from the file/new project option. In "Project Types" select "Visual C++ Projects." Under 
"templates" select "Win32 Project" and enter a name for the project (the name of the 
executable will be the project name with ".exe" appended). Click OK. Next the Win32 
Application Wizard opens. Click "application settings." Select "console application" and 
"empty project" then click on "finish." Now add all the source and header files to the project 
with the "project/add existing files" option. The debug and optimized versions are referred 
to as "debug" and "release." Compile the source with the "build" commands in the "Build" 
menu. 



4. Executing the Code 

If there were no errors in the previous step, you should now have a working executable. 
Next, simply type the executable filename at the command prompt in the directory where 
your exectuable was created. Note that some systems will not automatically look for the 
executable in the current directory. In this case, the command would be "./helium.exe". 
The code should respond by printing the following output and then exiting: 



USAGE: executable nmax Te topoff ST-mixing JobName 
executable - the executable filename 

nmax - the highest principal quantum number to explicitly consider 

Te - the electron temperature in Kelvin 

topoff - switch that enables (1) or disables (0) recombination above nmax 
ST-mixing - switch that enables (1) or disables (0) singlet-triplet mixing 
JobName - a name for this calculation, will also appear in output filenames 



-3- 



Each of the six items above must be specified and tfiey must be specified in that exact order. 
The minimum nmax is 15. Running the code with nmax > 100 is not recommended because 
the code is not well tested beyond that limit. (See Bauman et al. [2005] for a discussion of 
convergence with respect to nmax.) If singlet-triplet mixing is disabled, the results will be 
for the pure LS'-coupling case. Here is an example: 

. /helium. exe 15 10000 1 1 testrun 

This example will produce results for nmax=15 at a temperature of 10,000 K, with topoff 
and ^T-mixing included. The name "testrun" will be assigned to the job. 

5. Interpreting the Results 

The above example will produce the following output files: 



testrun. 


.15. 


.Output . txt 


testrun. 


.15. 


.Results .txt 


testrun. 


.15. 


.Results_Ryan . txt 


testrun. 


.15. 


.Levels_st .txt 


testrun. 


.15. 


.Thrshld_X_Scts . txt 


testrun. 


.15. 


.Solution_st . txt 


testrun. 


.15. 


.Vector.st .txt 


testrun. 


.15. 


.OscStrMethod_st . txt 



The first two parts of each output filename are given by the JobName and nmax parameters. 
The last part of the filename (before the ".txt") describes the file contents as shown in 
Table 1. 

6. Referencing Use of the Code 

The primary developer of the code was R. P. Bauman. Any publication that uses the 
results of the code should include a reference to Bauman et al. (2005). For reproducibility, 
the input parameters should also be reported. 
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Table 1. Output filenames and contents. 



Last Part of Filename 


Contents of File 


Output 


Basic sanity checks and various data; this file can be quite cryptic. 


Results 


J-resolved wavelengths, emissivitites, oscillator strengths, 




transition probabilities, occupation numbers, etc. 


Results-Ryan 


Same as Results but nLS'-resolvcd for comparison with 




the He I model in Cloudy (see Porter et al. 2005, 2007). 


Levels-st 


Various information about levels. The "st" in the filename 




will be "Is" if ST-mixing is disabled. 


Thrshld-X-Scts 


Threshold photoionization cross-sections from various sources. 


Solution-st 


The solved vector of level populations (in atomic units). 


Vectoi-st 


The initial vector of recombination rates (in atomic units). 


OscStrMethod-st 


The sources of oscillator strengths for each transition. 



